Cyclic Spin Exchange in Cuprate Spin Ladders 
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We investigate the influence of a cyclic spin exchange Jcyc on the one- and two-triplet excitations of 
an undoped two-leg 5=1/2 ladder, using the density matrix renormalization group (DMRG). The 
dispersion of the 5=0 two-triplet bound state is dramatically reduced by Jcyc due to a repulsion 
between triplets on neighboring rungs. In (La,Ca)i4Cu24 04i a consistent description of both the 
spin gap and the 5=0 bound state requires Jcyc/ J ± ~ 0.20-0.27 and J/Jx ~ 1.25-1.35. With these 
coupling ratios the recently developed dynamical DMRG yields an excellent description of the entire 
5=0 excitation spectrum observed in the optical conductivity, including the continuum contribution. 

PACS numbers; 75.10.Jm, 75.40. Gb, 75.40. Mg, 74.72.Jt, 75.30. Et 

The antiferromagnetic parent compounds of the high- 
Tc cuprates are thought to be the best representatives 
of the two-dimensional 5=1/2 square-lattice Heisenberg 
model. Understanding their magnetic properties is of ut- 
most importance due to the intimate relationship of mag- 
netic correlations and high- Tc superconductivity. Re- 
cently, the question how to set up a minimal model which 
accounts for these magnetic properties has been read- 
dressed. High-resolution inelastic neutron scattering ex- 
periments performed on the two-dimensional 5=1/2 anti- 
ferromagnet La2Cu04 0| exhibit a magnon dispersion at 
the zone boundary which cannot be obtained within a 
nearest-neighbor Heisenberg model. It has been argued 
that the inclusion of a cyclic spin exchange term of about 
20% would reproduce this dispersion [Q. This cyclic 
spin exchange emerges as a correction to the nearest- 
neighbor Heisenberg Hamiltonian in order t"^ /U^ from a 
i/L/-expansion of the one-band Hubbard model ||]. It is 
expected to be the dominant correction within a more 
realistic three-band description of the Cu02-planes be- 
cause there the cyclic permutation of 4 spins on apla- 
quette can take place without double occupancy 
Similar cyclic spin exchange processes have proven to be 
significant in other systems, e.g. in cuprate spin chains 
a ferromagnetic 2-spin cyclic exchange is responsible for 
the unusually strong exchange anisotropy 

Cuprate spin ladders offer an alternative approach to 
decide about the existence and potential implications of 
a cyclic spin exchange term. They are composed of the 
same corner-sharing Cu-0 plaquettes as the 2D cuprates, 
thus similar exchange couplings are expected for all spin 
products. In fact the inclusion of a cyclic spin exchange 
has been suggested in order to explain the smallness of 
the ladder spin gap observed in (La,Ca)i4Cu2404i 
However, it is impossible to extract a unique set of cou- 
pling parameters or even to settle the existence of a cyclic 
spin exchange term from an analysis of the spin gap only. 
Here, the optical conductivity cr(w) (l^Jlltl can provide 
the missing information. Magnetic excitations can be 
observed in (j{lo) via the simultaneous excitation of a 



phonon |12| , |13| . In (La,Ca)i4Cu2404i, two peaks in oiuj) 
were identified as 5=0 bound states of two triplets |10| . 

Here we show that a consistent evaluation of a{uj) is 
achieved by inclusion of Jcyc- The dispersion of the 5=0 
bound state is very sensitive to the magnitude of Jcyc- 
The combined analysis of both the 5=0 bound state and 
the spin gap establishes a unique set of exchange cou- 
plings, namely J^jc/Ji -0.20-0. 27, J/ w 1.25-1.35 and 
J_L = 950-1100 cm-i for (La,Ca)i4Cu2404i. For these 
values, we calculate the entire spectrum of <j{lo) using 
the dynamical DMRG |l^ and find excellent agreement 
with new experimental data which yield an improved es- 
timate of the continuum contribution. This clearly con- 
firms the above value for Jcyc- One might speculate that 
such a sizeable cyclic spin exchange term will have sig- 
nificant consequences for superconductivity in the doped 
spin ladders [ p^p^ since magnetic and pairing correla- 
tions are considered as closely related phenomena [ p^ . 

We study an isolated two-leg 5=1/2 ladder charac- 
terized by an antiferromagnetic Heisenberg Hamiltonian 
with an additional cyclic spin exchange term: 

H = J± ^ Si^iSi^r + J ^ (Si,;Si+i^; + Si_rSi+i,r) (1) 
i i 

+ Jcyc'^^(P(iJ)^(i^r),(i+l,r),(i+l,l)+P{i}){i,r){i+l,r){i+lA)^ 
i 

where J± and J denote the rung and leg couplings, the 
index i refers to the rungs, and I, r label the two legs. 
The cyclic permutation operator P1234 for 4 spins on a 
plaquette is given by: 



P 



1234 



p-1 _ 
-^^1234 ~ 



(2) 



4(SiS2)(S3S4) + 4(SiS4)(S2S3) - 4(SiS3)(S2S4) 
-I-S1S2 + S3S4 + S1S4 + S2S3 + S1S3 + S2S4 . 

The influence of Jcyc can easily be understood in terms 
of the representation of the spin operators by rung sin- 
glets s- and triplets 4q {a = x,y,z) p^jlgj ]: 



= x,y, 

2 Jcyc 



(3) 
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FIG. 1 . DMRG results for a N=80-site ladder with J = Jj_ 
and < Jcyc/ J < 0.3. (a) One-triplet dispersion. The inset 
shows the spin gap As (extrapolated to A'^ = oo j2^) as a 
function of Jcyc/ J together with an exponential fit. (b) Lower 
edge of the two-triplet continuum (open symbols) and the 
S = Q two-triplet bound state (full symbols). 
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The consideration of Jcyc primarily renormalizes the 
coupling strengths of the individual terms in the original 
Heisenberg Hamiltonian. The only new contribution is 
the last term, a repulsive interaction between triplets on 
neighboring rungs. With respect to the renormalization 
of the couphng constants the strongest effect is a reduc- 
tion of the triplet on-site energy. Consequently the spin 
gap as well as the energy of the entire one-triplet disper- 
sion are reduced. In Fig. |^(a) we show the one-triplet dis- 
persion, obtained by the Lanczos-vector method [ ^J2l| 
for a N=80-site ladder with J = J± and different values 
of Jcyc- The shape of the one-triplet dispersion remains 
qualitatively unchanged. The strongest suppression af- 
fects the minimum {px = tt) and the maximum (near 
Px W7r/3). Nevertheless, there still is a local minimum 
at Px—0 even for Jcyd J — ^-i, in contrast to the exact 
diagonalization results for a N=24-site ladder in Ref. [Q. 
Extrapolating to an infinite ladder j22j, we find that the 
spin gap = Eo{N^oo,S=l) - Ea{N = 00, S = 0) de- 
creases exponentially (inset of Fig. |l](a)), but we cannot 
conclusively decide from our present data, whether the 
spin gap indeed vanishes at Jcyc/ J ~ 1-2, as has been 
proposed in Ref. p4[ . 

Similar considerations also apply to the two-triplet ex- 



citations. With increasing Jcyc the lower edge of the 
two-triplet continuum (open symbols in Fig. 0(b)) shifts 
corresponding to the one-triplet dispersion. In addi- 
tion, however, a repulsive interaction between neighbor- 
ing triplets is introduced by Jcyc, as mentioned above (see 
Eq. (^)). This reduces the binding energy of the 5'=0 two- 
triplet bound state. As a consequence the width of the 
bound-state dispersion decreases dramatically (full sym- 
bols in Fig. |^(b)). Evidently the dispersion of the 5=0 
two-triplet bound state is very sensitive to the magnitude 
of Jcyc and thus provides a suitable quantity to probe the 
existence and magnitude of a cyclic spin exchange. 

In (La,Ca)i4Cu2404i the S=0 bound state has been 
observed in the optical conductivity cr{uj) via phonon- 
assisted two-triplet absorption |l^. Although the total 
momentum of the excitation has to be zero, the simulta- 
neous excitation of a phonon with arbitrary momentum 
permits to measure a weighted average of the S=0 exci- 
tations from the entire Brillouin zone. In Fig. |^ we show 
the magnetic contribution to (j^lo) in La5.2Ca8.8Cu2404i 
p5| . For polarization of the electrical field parallel to 
the legs, the two peaks in a\cg{u!) at u;i=2140 cm~^ and 
ti;2=2780 cm~^ correspond to van Hove singularities re- 
sulting from the maximum Wmax and minimum tJ,nin of 
the S=0 bound-state dispersion at Px~T^/'2 and Px = tt, 
respectively (see Fig. |l|(b)). The broad features at higher 
energies represent the continuum contribution. 

In order to determine the exchange couplings for 
La5.2Ca8.8Cu2404i we have calculated ujma.x and tUnun 
for different coupling ratios J/J± and Jcyc/ J l using 
the Lanczos-vector method [Q. By extrapolation to 
an infinite ladder we obtain the values displayed in 
Fig. ^(a). The corresponding results for the spin gap 
As are shown in Fig. ^(b). In the range depicted in 
Fig. ^ these values may be approximated within 5% by 
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FIG. 2. (a) Maximum of the S=0 bound-state dispersion 
at ~ 7r/2 (full symbols), minimum at = tt (open sym- 
bols) and (b) spin gap As as a function of J/Jx for different 
values of Jcyc/ J± ■ AH results have been extrapolated to N=cx3 
0. 
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!l.64J_L-0.54J, 



1.87Jcj^c+l-53J 



0.61J 

and As wO.48 J_L — 0. 84 Jcyc- In principle the three exper- 
imental quantities wi, loi and As ~ 280 cni^^ H] should 
suffice to determine J eye J and J ^. However, wi and loi 
depend also on the frequency Wp/i of the phonon involved 
in the optical absorption process, i.e., ijJ\=ijj^-^ 



Px—'^ 



and L02 = tjJmax + ^ph~^ ■ Making use of the relative 
insensitivity of As and c^min on J/J±, we can deter- 
mine the magnitude of the cyclic spin exchange with- 
out consideration of lu2 and hence omit any effect of a 
possible phonon dispersion. For a conservative estimate 
of the Cu-0 bond-stretching phonon frequency i^p^~^ = 
600 ± lOOcm-i we obtain Jcyc/J± = 0.20 - 0.27 with 
J_L w 950 — llOOcm"^. The dispersion of the phonon en- 
ters only in the determination of the coupling ratio J/J±. 



Allowing for a dispersion of 



t/2 



^ph - ±50 cm- 



we obtain J/J±_ ~ 1.25 — 1.35. 

Up to now we have derived a set of exchange couplings 
which consistently describes the spin gap and the S=0 
two-triplet bound state. Now we compare the spectral 
density calculated for this parameter set with the line 
shape of the entire spectrum observed in cr(u;). This 
provides a thorough test whether we indeed have found 
the minimal model which captures all relevant proper- 
ties. We consider an isolated Cu203-ladder, for which 
the major contribution to (j{uj) results from the simul- 
taneous excitation of two neighboring spin-flips plus a 
Cu-0 bond-stretching phonon: 



J2 fplm{{SB-p;SBp))^^_, 



(4) 
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FIG. 3. Symbols: Correction-vector results for the optical 
conductivity o-(a;) of an N=80-site ladder with polarization of 
the electrical field E along the legs (top panel) and along 
the rungs (bottom panel), using J/ J_l — 1.3, .leve l J ±=0.2, 
Jx^lOOO cm-\ w^^^570 cm-\ a;X'^=620 cm~^ H and a 



Px Py=Q,-IT 



finite broadening of 5 = O.lJx. The leg polarization contains 
two contributions, where the two legs are in-phase {py=Q) or 
out-of-phase [py = tt) with each other. Gray lines: (j(uj) of 
La5.2Ca8.8Cu2404i at T=4 K [^. Inset: one-triplet disper- 
sion (dashed line), lower edge of the two-triplet continuum 
(thin solid line) and S=0 bound state (thick solid line) for 
the above parameters. 



where <54<=*^=^E.E 



and (5B™"g 



{Si^iSi^r - {^i,i^i,r)) are the 



spin-flip operators for polarization of the electrical field 
along the legs and the rungs, respectively. To lowest or- 
der (4th order in the Cu-0 hopping tpd), the dominant 
contribution to the phonon form factor f}^'^ comes from 
the in-phase and the out-of-phase stretching modes of the 
0-ions on the legs, whereas for the out-of-phase 

stretching mode and the vibration of the 0-ion on the 
rung are taken into account: 



rung 



5sin (y) 



(5) 



So far we have used the Lanczos-vector method 
(LVM) to calculate the dispersion of the bound state. 
This was justified by the fact that the LVM uses the first 
Lanczos vectors as target states besides the ground state 
and consequently represents the bound state quite accu- 
rately. However, in order to calculate the high-frequency 
range of (t{uj) the LVM is not reliable. We use the 
correction- vector method [fi 26|, because here the den- 



sity matrix is optimized to represent the correction vector 
for each frequency point separately. 



In Fig. H we compare (j{w) for polarization along the 
legs and along the rungs, obtained by the correction vec- 
tor method for a ladder with N=80-sites, with the exper- 
imental result of La5.2Ca8.8Cu2404i |2^. Here we have 
used J/J^=1.3 and Jcyc/^±=0.2 with J± = 1000 cm~^ 
psf . For polarization along the legs, (t{uj) contains 
two contributions, where the two legs are excited in- 
phase {py = 0) or out-of-phase {py = tt). The in- 
phase mode contains two- and four-triplet excitations, 
as can be inferred from the rung-triplet representation of 
the spin-fiip operator Si^;Si+ij + Si^rSi+i^r in Eq. (|^). 
It includes the 5=0 two-triplet bound state as dis- 
cussed above. The out-of-phase mode reflects excita- 
tions of three triplets due to Si_/S,i+i_/ — S^.^S^+i^^ = 

and hence contributes only to the high-energy contin- 
uum excitations |p9|| . For polarization along the rungs, 
cr{uj) contains only two-triplet excitations (see Eq. (||) for 
the rung-triplet representation of SijSi^r) and the lower 
bound state is suppressed due to a selection rule . 

Comparing the DMRG results for polarization along 
the legs with the experimental spectra we find excellent 
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agreement within the range of the bound state (which 
has been used in order to fix the coupling constants), 
and also with respect to the overall line shape of the 
high-energy continuum excitations. For this polariza- 
tion the experimental determination of the precise spec- 
tral weight of the continuum is difficult due to an elec- 
tronic background | p0| , pl]j25| . The small differences be- 
tween experiment and theory in the continuum range are 
certainly within the experimental error bars. For polar- 
ization along the rungs, where experimentally the mag- 
netic contribution to (t{uj) can be identified unambigu- 
ously and with high precision |jl^,|ll|, the continuum is 
reproduced nearly perfectly. This provides strong evi- 
dence that a nearest-neighbor Heisenberg Hamiltonian 
with J/J±=1.3 and an additional cyclic spin exchange of 
Jcyc/ J±=0-'2' represents the magnetic excitations of the 
spin-ladder compound (La,Ca)i4Cu2404i. 

In spite of the overall impressive consistency of exper- 
iment and theory, some discrepancies are visible: (i) In 
CTicg(w) the upper bound state (for « ^/2) is sup- 
pressed much stronger in the DMRG calculation than in 
the experimental spectrum, (ii) In crrung('^) the exper- 
imental data displays a pronounced shoulder at about 
2200 cm~^, i.e., shghtly above the suppressed lower 
bound state, which is not present in our DMRG calcu- 
lation. We propose that these discrepancies are related 
to different phonon form factors which arise due to the 
actual arrangement of the ladders within the trellis lat- 
tice, and in higher order (with respect to Cu-0 hopping) 
in the external photon-electron-phonon vertex. Consid- 
ering, e.g., an additional term oc sin^(pa;/2) in will 
enhance the weight of the upper bound state in aicg{uj) 
with respect to the lower one. Such considerations, how- 
ever, depend on details of the non-local excitation process 
in these ladder systems and are therefore not of funda- 
mental relevance. 

In summary, we have demonstrated that the inclu- 
sion of a cyclic spin exchange Jcyc is necessary to ob- 
tain a consistent interpretation of optical experiments 
and the magnitude of the spin gap in (La,Ca)i4Cu2404i. 
The dispersion of the 5*^0 bound state is very sensi- 
tive to the magnitude of Jcyc- Quantitative investiga- 
tion yields Jcyc/J± ~ 0.20-0.27 and J/J± « 1.25-1.35 
with J± « 950-1100 cm~^. With these parameters we 
obtain excellent agreement with the experimental spec- 
tra of La5.2Ca8.8Cu2404i, both for the bound state and 
the continuum. This gives strong evidence that we in- 
deed have identified the minimal model which captures 
all relevant magnetic properties. 
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